Caging dynamics in a granular fluid 



P.M. Reis*, R.A. Ingale and M.D. Shattuck 

Benjamin Levich Institute, The City College of the City University of New York 
UOth St. and Convent Av., New York NY 10031, USA 

We report an experimental investigation of the caging motion in a uniformly heated granular 
fluid, for a wide range of filling fractions, <j). At low 4> the classic diffusive behavior of a fluid is 
observed. However, as <j) is increased, temporary cages develop and particles become increasingly 
trapped by their neighbors. We statistically analyze particle trajectories and observe a number of 
robust features typically associated with dense molecular liquids and colloids. Even though our 
monodisperse and quasi- 2D system is known to not exhibit a glass transition, we still observe many 
of the precursors usually associated with glassy dynamics. We speculate that this is due to a process 
of structural arrest provided, in our case, by the presence of crystallization. 
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In additional to being of great industrial and geological 
importance, granular materials are of fundamental inter- 
est due to their strong non-equilibrium nature En- 
sembles of granular particles are intrinsically dissipative 
and any dynamical study must involve energy injection. 
These ingredients make the understanding of granular 
assemblies a challenging endeavor and a general theoret- 
ical framework is still lacking. A possible approach is the 
study of non-equilibrium steady states, in systems where 
energy injection perfectly balances dissipation. Whereas 
some progress has been made in the fast dilute regime 
0, the understanding of the dense case of these steady 
states f remains an open question. One avenue of research 
has been to borrow concepts from other dense assemblies 
of particles such as colloids, suspensions and emulsions. 
Moreover, it is thought that connections with molecular 
glass formers may be relevant |2j . A priori, it is far from 
obvious that such analogy may be valid due to the enor- 
mous differences in lengthscales and the the mechanism 
of energy supply. However, one common feature in all 
of these seemingly disparate systems is the presence of 
cages: each particle is temporarily trapped by its neigh- 
bors and then moves in short bursts due to nearby coop- 
erative motion. This results in highly heterogeneous mo- 
tion and slowing down of the dynamics. In molecular sys- 
tems this behavior is typically observed indirectly from 
scattering experiments [4]. In colloids, however, caging 
motion has been observed directly through microscopy, in 
both 3D and quasi-2D Q geometries. A large number 
of theoretical 0, E| and numerical 0, [H3 studies have set 
out to further investigate the importance of this hetero- 
geneous dynamics. The relevance of the caging in driven 
granular materials 0, and air-fluidized particle sys- 
tems 43] has only recently started to be addressed. In 



particular, Dauchot et. al. [3 have reported on a gran- 
ular system driven by cyclic shear where they observed 
many "glasslike" features, but for a single value of the 
filling fraction. 

In this Letter we present a novel dynamical study of 
a monodisperse, quasi-2D granular fluid Q, in which 
the particles are excited by a spacially uniform stocastic 
forcing. In our system the filling fraction, </>, can be var- 
ied from a single particle to hexagonal close packing. We 
have shown that the structural configurations of this non- 
equilibrium stea dy s tate are the same as those of equilib- 
rium hard-disks [l4[. In particular, the system exhibits 
three phases: an isotropic fluid phase <\> < <\>\ = 0.652, 
a crystaline solid phase <\> > (j) s = 0.719, and an inter- 
mediate^ phase (j>i < <j) < (j) s consistent with a hexatic 
phase [l^. The phase boundaries, <pi (liquidus point) and 
4> s (solidus point) are determined by structure alone Q. 
Here, we explore the dynamics of these phases through 
single particle trajectories, focusing on the caging dy- 
namics seen in the intermediate phase. In Fig. we 
present typical single particle trajectories for filling frac- 
tions in each of the three phases. Simple fluid behav- 
ior is observed at low 0, characterized by random diffu- 
sion (Fig. HJl). Above crystallization (</> > <j) 3 ) particles 
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FIG. 1: Experimental frames with superposed typical trajec- 
tories of a single particle: (a) <j) — 0.567, (b) <j) — 0.701 and 
(c) <j> = 0.749. Note that even though only a single trajec- 
tory is shown for each 0, particle tracking and statistics were 
performed over all particles within the imaging window. The 
scalebar is 2mm. 
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become fully arrested by thier six hexagonally packed 
neighbors (Fig. [IJ). In the intermediate phase, we see a 
mixture of these behaviors (Fig. At short times, par- 
ticles are temporarily trapped in cages formed by thier 
neighbors, but at long times they diffuse from cage to 
cage. We will use the Mean Square Displacement (MSD) 
and the Intermediate Scattering Function (ISF) to show 
that the caging dynamics seen here is qualitatively identi- 
cal to that of dense molecular and colloidal systems 0, Q 
and supercooled liquids [H| . This is surprising since our 
experiment is fundamentally different. The associated 
lengthscales of our granular system are typically 3-7 or- 
ders of magnitude larger than those of molecular systems 
and colloids. Moreover, the steady states we study are 
inherently far from equilibrium since energy is both in- 
jected (through vibration) and dissipated in inter-particle 
inelastic collisions and frictional contacts. 

Our experimental apparatus consists of vertically vi- 
brated, D = 1.19mm diameter stainless steel spheres 
confined between two 85. 3D diameter horizontal glass 
plates, separtated by 1.6D, which constrain the parti- 
cle motion to be quasi-2D. Our system is described in 
more detail in UA UM and improves on Olafsen and Ur- 
bach's similar system by using a roughned bottom 
plate. This allows us to study a wide range of filling 
fractions (1.4 x 10 -4 < (j) < 0.8). We sinusoidally vi- 
brate the system with a frequency / = 50Hz and a max- 
imum acceleration equal to 4 times gravity, but the struc- 
ture and dynamics of the system are independent of the 
forcing for a wide range of parameters [l^. To ensure 
repeatable initial conditions, we start with all particles 
hexagonally packed near the boundary. A steady-state is 
reached by waiting for 12000 vibration cycles before 10 
seconds of data are acquired. We record the dynamics in 
a (15 x 15mm 2 ) central region using a high-speed camera 
at 840 Hz and track the trajectories of all particles. 

We first analyse the particle trajectories by measuring 
the Mean Square Displacement defined as, 

M(t) = ([r(t)-r(0)} 2 ), (1) 

where r(t) is the position of a particle at time t, r(0) 
is its initial position, the brackets () signify ensemble 
averaging over many realizations and time invariance is 
assumed. The MSD for a range of <p are shown in Fig. 
El For the case of a single particle in the cell (marked 
with a square in Fig. [2j, the motion at short times is 
ballistic, and M(t) ~ t a where a ~ 2 (a = 2 for pure 
ballistic motion). At later times, the particle moves dif- 
fusively and the slope of M{t) tends to a = 1. This 
shows that the trajectory of a single particle is indeed 
randomized across the cell. For all (/>, the motion at early 
times is superdiffusive with a ~ 2 showing ballistic mo- 
tion. For <j) < 0.719 the motion always becomes diffusive 
at long times with M(t) ~ t. Eventually, this becomes 
increasingly noisy due to the lack of statistics in the time 
averaging. For <j) > 0.719 the particles are trapped by 
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FIG. 2: Time dependence of the Mean Square Displacement 
for various values of filling fraction (numerical values given in 
the box) . The curve marked with a square is for a single parti- 
cle in the cell. The arrow points in the direction of increasing 
(j). Along the arrow, the symbols (*) and (+) are located at <j)i 
and S , respectively. The horizontal line corresponds to the 
area of l/4th of the system, above which finite system size 
effects become important. 



their six hexagonally packed neighbors, and M(t) levels 
off to a constant value set by the lattice spacing. In the 
intermediate phase, however, a plateau emerges at in- 
termediate times where the motion is subdiffusive with 
< a < 1. This plateau becomes obvious above the liq- 
uidus point, <pi (marked as * in Fig. [3J) and represents 
the slowing down due to the cage effect shown in Fig. 
A similar dependence of M(t) has recently been observed 
in a quasi-2D system of bidispersed particles fluidized by 
a uniform upflow or air |l3|. 

Another classic measure in the study of dense liquid 
phases is the Intermediate Scattering Function Q which 
is defined as, 

^(q,*) = j; I> x p Hq • [rj(t) - Tj(p)])), (2) 
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where q is a wavenumber and Tj(t) is the trajectory of 
particle j out of TV particles in the system. This mea- 
sure is widely used in colloids since it is readily available 
through light scattering experiments 4] and is, essen- 
tially, a measure of the time decorrelation of the posi- 
tional wavevectors. In dense colloids and supercooled 
liquids, F s (t,q) captures the relaxation due to caging in 
the form of a two-step relaxation: 1) the fast (early time) 
(3 relaxation which corresponds to the diffusion inside the 
cage followed by 2) the a relaxation corresponding to the 
time it takes for the particle to diffuse out of the cage. 

In Fig. 01 we plot F s (t) for wavevector qD = 2.14 
(q = 1.8mm" 1 ) for various values of (j). As expected, in 
the crystal phase (<p > 0.719) F s levels at a value close 
to 1 and little decorrelation occurs, since each particle is 
fully trapped. In the fluid phase, F s rapidly decays as 
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FIG. 3: Time dependence of the Intermediate Scattering 
Function, with qD — 2.14, for various filling fractions (nu- 
merical values of <j) given in the box). The arrow points in 
the direction of increasing <j>. Along the arrow, the symbols 
(*) and (+) are located at <j>i and <j> s , respectively. The solid 
lines are fits to Eqn. @. 



particles move across the cell diffusively and the initial 
positional wavevector quickly decorrelates. In the inter- 
mediate phase we observe the classic a and /3 two-step 
relaxation; there is a clear intermediate plateau and the 
a relaxation occurs at increasingly longer timescales, as 
cj) is increased. As for the MSD before, this two-step 
relaxation becomes visible above the liquidus point, <\>\ 
(marked as * in Fig. |3J). For each value of (/>, the a 
relaxation at late times is well described by a stretched 
exponential of the form, 



F s (q,t)~eM- (t/r(q)f iq) ] 



(3) 



where r(q) is a relaxation time and the stretching expo- 
nent is typically f3(q) < 1. Fits of this stretched expo- 
nential form to the experimental data are shown as solid 
lines in Fig. [3J It is interesting to note that this behavior 
is in good agreement with the predictions of Mode Cou- 
pling Theory (MCT) 7]. This is highly surprising since 
MCT has been developed for thermal fluids and is not 
known to apply to non-equilibrium systems such as ours. 

We now focus on the (/-dependence of both the relax- 
ation time r(q) (Fig. Et) and the exponent /3(q) (Fig. 
Eb), for various values of <p. We defined r(q) as the time 
it takes for the experimental curves of F s to fall to a 
level of 1/e, i.e. F s (r) = 1/e. The stretching exponent 
(3(q) is the local slope of the quantity log(— log(F s )), in 
the neighborhood of r. At low filling fractions, r scales 
as q~ 2 . As the filling fraction is increased this scaling 
continues to hold but only up to a cut off value above 
which (lengthscale below which) it sharply drops. This 
cutoff lengthscale can be associated with a characteristic 
size of the cage, which becomes increasingly smaller as <fi 
is increased. On the other hand /?, the local exponent, 
tends to one at small q but decreases progressively below 



one for higher filling fractions and F s becomes increas- 
ingly stretched. These findings can now be combined 
and interpreted as follows. For small q (i.e. large length- 
scales), r ~ q~ 2 and /3(q) — > 1, which together imply 
that F 3 (q,i) ~ exp(— Dq 2 t). This is the result for clas- 
sical diffusion 0. For large q (i.e. small lengthscales) 
this Brownian scaling breaks down to a stretched expo- 
nential with /3 < 1, which can attributed to the presence 
of dynamic heterogeneities due to caging. 

Returning to the case of fixed qD = 2.14, we plot r 
as a function of (j) in Fig. A significant slowing down 
of the dynamics can be seen at high 0, as crystallization 
is approached. It is highly surprising that this slowing 
down with <p is well described by the Vogel-Fulcher law 
[20| as found in many glass forming systems, 



exp[A/(</> c 



(4) 



where A = 0.094 ± 0.004 is a fitting parameter but 
<j) c = 0.719 ±0.007 is the filling fraction for crystallization 
which was determined independently from experiments 
Q. A power-law fit was not as satisfactory. Note that, 
as the fluid goes through the transition from isotropic 
fluid to the intermediate phase (at <\>\ = 0.652), r shows 
no particular feature. This functional dependence of the 




FIG. 4: (a) Wavevector dependence of the relaxation time, 
r and (b) local stretching exponent /3, for various values of 
filling fraction. The arrows point in the direction of increasing 
(j) and the numerical values of (j) are given in the boxes. Along 
the arrow, the symbols (*) and (+) are located at (j>i and <j> s , 
respectively. 
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FIG. 5: Relaxation time, r, extracted from the Intermediate 
Scattering Function as a function of filling fraction. The solid 
line is a fit to the Vogel-Fulcher law of Eqn. (|p). Dashed and 
solid lines represent the location of the liquidus and solidus 
points, respectively. 



relaxation has also recently been found within the gran- 
ular context by Fierro et. at. 21] in a numerical lattice 
model. 

In summary, we have studied the dynamics of a uni- 
formly heated granular fluid. We have observed a number 
of features typically associated with dense liquid behavior 
in molecular and colloidal systems, namely: prominence 
of cages, development of a plateau in the MSD and ISF 
with the breakdown of the Brownian diffusive behavior 
and a Vogel-Fulcher relaxation. These features become 
particularly visible above the liquidus point, for <j)> <j>i. 
In particular, our results can be directly compared to 
both experiments _6] and simulations ^(j of quasi-2D 
colloidal particles. This is surprising since our exper- 
iment is fundamentally different: it is intrinsically far 
from equilibrium since energy is not conserved, the con- 
stituent particles are macroscopic and it is known that for 
monodisperse and quasi-2D systems such as ours there is 
no ideal glass transition (the phase diagram is fluid, 
intermediate phase and crystal). In the absence of a 
structural glass, we propose that the spatially uniform 
stochastic way of injecting energy along with a process 
of structural arrest provided by crystallization |T| are 
at play to account for the many observed similarities. 
This suggests that theoretical frameworks previously de- 
veloped for dense thermal liquids, for example MCT Q, 
might shed some light to the description of excited gran- 
ular materials. 
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